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ABSTRACT 

Using HST images, we separate the bulge-like (dubbed pbulge) and disk-like (pdisk) components of 71 
galaxies in the rich cluster MS 1054-03 and of 21 in the field. Our key finding is that luminous pbulges are very 
red with restframe U -B ^QA5, while predicted colors are bluer by 0.20 mag. Moreover, these very red colors 
appear to be independent of environment, pbulge luminosity, pdisk color, and pbulge fraction. These results 
challenge any models of hierarchical galaxy formation that predict the colors of distant {z ^ 0.8) luminous field 
and cluster bulges would differ Our findings also disagree with other claims that 30% to 50% of bright bulges 
and ellipticals at z ~ 1 are very blue {U-B < 0). 

Subject headings: cosmology:observations - galaxies:photometry - galaxies: fundamental parameters - galax- 
ies:evolution - galaxies:formation 



1. INTRODUCTION 

The ages of bulges (also commonly referred to as 
spheroids)^ remain an important unsolved problem in stellar 
populations and galaxy formation (Wyse, Gilmore, & Franx 
1997). The three current and contending paradigms are early, 
rapid single-collapse formation of bulges and ellipticals at 
redshifts z > 2; formation at later epochs via merging of disk 
galaxies; and the formation of bulges through secular evolu- 
tion processes in disks (see review by Kormendy & Kennicutt 
2004). In general, distant cluster galaxies show consistent re- 
sults from different studies, namely a tight color-magnitude 
relation for E-SO and very red colors that persist to redshifts 
z '--^ 1 (see, e.g., van Dokkum et al. 2000, Blakeslee et al. 2003, 
Menanteau et al. 2004). Both support the early bulge forma- 
tion scenario. In contrast, field** E-SO and bulges show incon- 
sistent and inconclusive results. For example, using the lumi- 
nosity function as a test between the two scenarios, Lilly et al. 
( 1 996) find little luminosity or density evolution for red galax- 
ies; Im et al. (2002 [GSSIO]) see significant luminosity evo- 
lution and little volume density evolution of morphologically 
selected E-SO to z ~ 1 ; and Bell et al. (2004) find luminosity 
and density evolution of red galaxies (presumed to be dom- 
inated by E-SO galaxies) that maintain a constant luminosity 
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density with lookback time to z 1 (see recent review by Mc- 
Carthy 2004). Perhaps the picture is mixed for galaxies with 
different masses, as suggested by fundamental plane studies at 
these redshifts showing less massive galaxies appearing to be 
younger (e.g., Treu et al. 2005). Using colors as a surrogate 
for estimated ages, some studies support the merger picture 
by finding or claiming a large fraction (30% to 50%) of field 
bulges and/or early-type galaxies at intermediate redshifts to 
be blue or to show evidence of recent star formation (e.g., Ko- 
dama. Bower, & Bell 1999; Abraham et al. 1999; Schade et 
al. 1999; Tamura et al. 2000; Menanteau et al. 2001 & 2004; 
Ellis, Abraham, & Dickinson 2001), while a few other sur- 
veys find uniformly red E-SO (GSSIO) or bulges (Koo et al. 
2005 [GSS8]) at the same epochs. 

As shown by Menanteau et al. (2001, 2004), the properties 
of field bulges versus cluster bulges provide powerful diag- 
nostics of formation mechanisms. They found, contrary to 
expectations of hierarchical formation models, that the bulges 
of field spirals appeared bluer on average than field ellipti- 
cals (c.f., GSS8 where both are found to be the same). Using 
cluster ellipticals as a benchmark for old stellar populations, 
Menanteau et al. (2001) estimated that "at z '--^ 1 about half the 
field spheroidals must be undergoing recent episodes of star 
formation." Moreover, as emphasized by the work of Benson 
et al. (2002), the resolved colors of spheroids are particularly 
effective probes of their formation history. 

Based on decomposing galaxies into bulge-like and disk- 
like components. Figure 5 in GSS8 shows that the luminous 
field pbulges are redder and have a shallower color-magnitude 
relation (CMR) slope than that derived from integrated colors 
of early-type galaxies in the MS 1054 cluster. Without remov- 
ing the contamination by potentially bluer disks, we do not 
yet know whether cluster bulges are indeed bluer than field or 
have a steeper CMR slope. Moreover, the bulge colors of the 
cluster spirals, i.e., non-early-type, have yet to be measured. 
Photometry that separates bulges from disks of cluster galax- 
ies is clearly needed. This letter presents the results of such a 
study. 

We adopt a Hubble constant Ho = 70 km s"'Mpc"' and a 
flat cosmology with = 0.3 and JIa = 0.7. Our photometry 
is in the Vega system with V referring to Feoe and / to f8i4 for 
the WFPC2 data from HST. At redshift z 0.83, V-I is very 
close to restframe U-B and / is close to restframe B, thus 
making any K-corrections to these measures very small. Ob- 
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served colors ofV—I ^ 2.2 are close torestframe U—B ^ 0.45 
and, for such very red stellar populations, / ~ 22 is equivalent 
to blue absolute magnitudes of Mb ^ -20.56. 

2. OBSERVATIONS 

Cluster MS 1054-03 at redshift z ^ 0.83 was chosen for 
analysis since it is among the best observed distant clusters. 
The cluster has relatively deep (6600s in each of V and /) 
and extensive (6 pointings of WFPC2 covering about 30 □' 
of the region) HST imaging from HST program GO 7372 
(see Hoekstra et al. 2000, van Dokkum et al. 2000). It also 
possesses over 100 spectroscopic redshifts (Tran et al. 2005, 
in preparation). Flat fielding and cosmic-ray removal of the 
HST images are described in detail by van Dokkum et al. 
(2001) who also measured colors and magnitudes via aper- 
ture photometry for 8 1 galaxies with spectroscopic redshifts 
(van Dokkum et al. 2000). 

Since our goal is to isolate the bulge component, the re- 
duced HST V and / images were analyzed with our own pho- 
tometry package (GIM2D: see Simard et al. 2002[GSS2] and 
GSS8 for details) to separate the luminosities and colors of the 
bulges and disks. GIM2D fits each object in an image with a 
PSF-convolved 2D two-component (bulge-ndisk) model. The 
two components are 1) the bulge, which is assumed to fol- 
low an r'/'* profile (de Vaucouleurs profile) that characterizes 
many luminous ellipticals and bulges (Andredakis, Peletier, 
& Balcells 1995), but this is a simplification (see de Jong 
et al. 2004 for a more in-depth study of local cluster bulge 
profile shapes using GIM2D), and 2) the disk, which is fit 
with an exponential profile. Since we cannot be sure that 
the components extracted from the decomposition are gen- 
uine bulges and disks, we referred to these measured compo- 
nents as photometrically-derived bulges and disks or photob- 
ulge and photodisk in GSS8. In this Letter, we dub these as 
pbulge and pdisk. 

Previous work with GIM2D (GSS2, GSS8) have shown that 
reliable two-band colors for compact bulges are derived by 
simultaneous fits of the profile parameters and the adoption 
of the simplifying assumption that color gradients for either 
of the subcomponents are small. One image in each filter is 
used in GIM2D. The pbulge and pdisk centers, sizes, orien- 
tation, and shapes are constrained to be the same for the two 
filters, while the pbulge-to-total ratio and total luminosity are 
left to float. Undersampled WFPC2 images with sub-pixel 
shifts cannot be re-sampled and stacked without loss of infor- 
mation, due to aliasing, unless a large number of images are 
available for stacking. Given that this is not case for our data, 
we make use, instead, of GIM2D's option of simultaneously 
fitting a series of non-rotated images in a single filter, each 
with its own set of PSF's. However, we cannot yet undertake 
simultaneous fits with multiple images in multiple filters. 

The cluster has 3 pairs of offset images in each of V and /. 
As a workaround to the limitations of GIM2D, we first fit each 
object with the 3 pairs of / images. The / was chosen since 
the bulges were generally very red and thus had higher S/N 
in / than in V. Except for the luminosity and pbulge-to-total 
ratio, all the other measured parameters from / were locked 
when the V image was fit. The conversion of the measured 
/ and y - / to restframe Mg and U-B were made using the 
method described by Weiner et al. (2005: GSS3) and Willmer 
et al. (2005), but yield values very close to that adopted by 
van Dokkum et al. (2000). Color conversion errors can be 
seen to range from 0.05 to 0.08 mag in Fig. 12 of Gebhardt et 
al. (2003: GSS9). 



The field bulges were chosen from a sample of galaxies 
within the V and / HST images of the "Groth Strip Sur- 
vey" (GSS), a "chevron strip" of 28 WFPC2 pointings ori- 
ented NE to SW at roughly 14:17h-52 at Galactic latitude 
b ^ 60deg (programs GTO 5090, 5109) covering an area of 
127 arcmin^ (Vogt et al. 2005). The galaxies were constrained 
to have about the same redshifts as MS 1054-03, namely z be- 
tween 0.79 and 0.87, where the spectroscopic redshifts were 
measured by DEEP'^. From a total sample of over 600 red- 
shifts (Weiner et al. 2005), 76 fell within our redshift range 
and had HST imaging. Except for one ultra-deep pointing 
(14:17.5h-52.5, Westphal Deep Survey Field 2) with 24,400s 
in V and 25,200s in /, the other HST images in GSS were 
comprised of a combined image in / of 4400s and in V of 
2800s. As detailed in GSS2, GIM2D processing of the 2800s 
V and 4400s / were made separately and simultaneously. To 
match the procedure for the cluster photometry using GIM2D, 
we adopted and locked the parameters from the 4400s / (ex- 
cept for luminosity and pB/T ratio) and used these to process 
the multiple images in V. The conversion of the input / and 
resultant V -I to restframe Mb and U-B were as mentioned 
above. 

Our prior experience with GIM2D (GSS8) suggests four 
additional criteria to define the final sample, with brackets en- 
closing the number of galaxies excluded by each criterion for 
the field and cluster, respectively: 1) [45,33] the pbulge had 
to be brighter than I = 23 (Mg ^ -19.5) to yield reliable pho- 
tometry; 2) [1,2] the pbulge half-light sizes were larger than 
one third of a pixel, i.e. 0.03 arcsec, again to assure reliable 
photometry; 3) [5,4] the pbulges had to have half-light sizes 
smaller than the half-light sizes of the pdisk subcomponent 
or the pbulges were dominant with more than 2/3 of the to- 
tal luminosity, i.e., pBAT > 0.67; and 4) [1,4] the morphology 
was not so distorted as to be poorly represented by our as- 
sumed two components. The third criterion was needed to 
reject cases where GIM2D fits a central probable bulge with 
a small exponential component and the outer probable disk 
with an r'/** component. Such structures appear to exist in na- 
ture but we chose to reject such anomalous bulges from this 
survey. Examples of rejected galaxies based on morphology 
include, e.g., a cluster galaxy that appears to be a spiral with 
two holes (ID 1403: see image in van Dokkum et al. 2000) 
and a field galaxy that appears to be a highly distorted close 
merger (GSS ID 073_1809: see image in GSS8). Another that 
was retained looks like a tulip with a stem (GSS ID 174_4356: 
see GSS8) and should perhaps be rejected; it is the lone very 
blue outlier in Fig. 1. Seven other blue pbulges seen in Fig. 
5 of GSS8 (most are not genuine bulges or have poor fits) lie 
outside the redshift range adopted here. 

Moreover, three of the field pbulges from GSS8 did not 
converge with our revised GIM2D procedures (GSS IDs: 
143_2770, 152_5050, and 292_6262). After the resti'ictions, 
the full sample of 1 14 MS 1054 cluster members was reduced 
to a final sample of 71 and the full GSS field sample of 76 was 
reduced to 21 pbulges. 

3. RESULTS 

The color magnitude diagrams (CMDs) of the pbulges are 
displayed in Fig. 1. To equalize the samples, the cluster 
pbulges have been segregated by their distance from the clus- 
ter center at J2000 10:56:59.93 -03:37:36.6 of galaxy 1484 in 

^ D eep Extragalactic Evolutionary Probe: see URL 

|http://deep .ucolick.org/ 
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van Dokkum et al. (2000). The virial radius or /?2oo measured 
as 241" by Tran et al. (1999) or 1.83 Mpc has been revised us- 
ing a larger sample of 130 galaxies and is 216" or 1.64 Mpc. 
All the cluster galaxies in our sample lie well within Rioo- 

The main result is clear from Fig.l: the slope, color, and 
dispersion of the CMDs are similar among the four sam- 
ples. Relative to the solid line showing the median color of 
U-B = 0.438 from the 71 cluster pbulges, the field pbulges 
appear slightly redder. We also see from Fig. 1 that the pbulge 
colors are red, independent of the distance from the cluster 
center (i.e., average density), pbulge to total (pB/T) ratio, and 
the color of the pdisk. We also find no support for claims or 
findings (see references in introduction) of 30% to 50% frac- 
tions of bulges or spheroidals being very blue (C/-B < 0) in 
the field at this redshift or within a broader range of redshifts 
of 0.7 to 1.0 (GSS8). Finally, we find that almost all pbulges 
are as red or redder than pdisks in both field and cluster. 

Before continuing, a few caveats and possible weaknesses 
of the survey are noted. The results here apply only to 1) 
bulges that are luminous (Mb ^ -19.5) and thus opens the 
possibility that lower luminosity bulges may have experienced 
a different formation history; 2) galaxies that are dominatd by 
one or two components, with the bulge having an r'/^ pro- 
file while the disk has an exponential; 3) galaxies that do not 
have large enough color gradients in either subcomponent to 
change the results; and 4) bulges that do not belong to obvi- 
ously merging systems. Our field sample with only 21 galax- 
ies is small, so that different properties might apply to 20% 
minorities. 

To quantify the CMD, Table 1 gives the results of fits to 
the data shown in Fig. 1 , usin g the biweight method of 
Heers, Flvnn, & GebhardtJ (Il990h . Since we found no lumi- 
nosity dependence of the colors, we locked the slope to and 
measured only the colors and dispersions. We confirm what 
is seen by eye, namely, that the cluster and field pbulges show 
very similar colors, with a hint that the inner cluster pbulges 
may be redder than the middle or outer cluster pbulges and 
that field pbulges may be slightly redder than cluster pbulges 
on average, but none of these differences are statistically sig- 
nificant (95% confidence limit or better, i.e., twice the RMS 
from Table 1). The color dispersions are also similar. Since 
the measured dispersions from the fits are smaller than that 
expected from the GIM2D color measurement errors, no reli- 
able intrinsic dispersions can be estimated, though the impli- 
cation is that they must be small. We estimate intrinsic color 
dispersions of less than 0.05 mag based on the variation of the 
measured values and GSS8 results. 

4. DISCUSSION & SUMMARY 

The main result is that the restframe colors, slope, and color 
dispersion of luminous (Mb < -19.5) cluster and field pbulges 
at redshift z ^ 0.8 are nearly the same. Making the simplify- 
ing assumption that colors are a surrogate for ages, we thus 
find that field and cluster pbulges are universally old (> 1.5 
Gyr) with no evidence that the ages of luminous pbulges 
are dependent on radial distance from the cluster center or 
whether they reside in the cluster or in the field, i.e. environ- 
ment. Moreover, while we do find slightly higher proportions 
of galaxies with large pbulge fractions (pB/T > 0.4) and with 
red pdisks (U-B > 0.25) in cluster MS1054 than in the field, 
we find that the very red pbulge colors in both the field and 
cluster are independent of pbulge luminosity, pB/T, as well as 
pdisk colors. 

The universally red pbulge colors with small dispersion 



found here can be contrasted with prior studies (see references 
in introduction) claiming that large fractions (30%-50%) of 
distant early-type galaxies have blue colors or large internal 
color dispersions that suggest active or recent star formation. 
These other studies differ from that undertaken here, either 
by including the underlying disk light in the photometry of 
the bulge or by including lower luminosity bulges. Our re- 
sults also disagree with models that predict cluster bulges are 
older than field bulges (e.g., Kauffmann 1996) or that bulge 
ages depend on bulge fraction (e.g., Kauffmann 1996) or on 
cluster-centric distance (e.g., Diaferio et al. 2001). 

Our findings for distant gala xies are, however, f ully consis- 
tent with the conclusions of P eletier et al. I (Il999h for a local 
sample of bulges in early-type galaxies. They derived from 
HST WFPC2 and NICMOS images the half-light colors of 
bulges in groups and the field. The colors were similar to that 
of elliptical galaxies in the Coma cluster and implied an age 
of 10 Gyr, corresponding to redshifts z ^ 2. The color spread 
within an effective radius where dust extinction was negligi- 
ble was also found to be equally small (^ 0.10 mag in B-I, 
equivalent to 0.05 mag in U-B), implying an age spread of at 
most 2 Gyr. Finally, because the bulge ages were so similar 
and old, they concluded that secular evolution of disks would 
have difficulty in forming bulges in early-type galaxies. Ex- 
ponential bulges of many late-type spirals, however, appear 
to have colors and thus presumably ages that are correlated 
with their disk colors, and lending support to secular evolu- 
tion scenarios for this class of bulges (Carollo et al. 2001; 
Kormendy & Kennicutt 2004). This work places no constraint 
on such bulges if they are of luminosity fainter than our limit 
of Mb - -19.5. 

The second key result is that the pbulge colors are very red 
with U-B - 0.45. This color lies between the U-B = 0.52 
at Mb = -20.5 of a fit to the CMD of 409 E-SO's from the 
RC3 (Schweizer & Seitzer 1992) and the U-B = 0.42 of 379 
SO's (Table 2 of Fukugita et al 1995) or U-B = 0.40 of 30 
bulges in early-type spirals (Peletier & Balcells 1996). The 
rough constancy of colors over the last 7 Gyr is a surprise 
for which we have no simple, compelling explanation if only 
single-burst stellar populations with standard IMF's are con- 
sidered. Since bulges should be bluer by 0.2 mag or more in 
U-B at lookback times corresponding to redshifts z ^ 0.83, 
we expect to observe typical colors nearer U-B ^ 0.20 - 0.32, 
depending on the formation epoch of bulges. If, however, an 
old, metal-rich stellar population were to be combined with 
small amounts of additional star-formation (or equally old but 
bluer stars of lower metallicity), significant changes in lumi- 
nosity from passive evolution of the stellar population can 
be accompanied by nearly constant restframe U-B colors 
(GSS8, GSS9; Schiavon et al. in preparation; Barker et al. 
in preparation). Such "drizzling" or "frosting" pictures find 
support ranging from detailed studies of local old stellar pop- 
ulations (Trager et al. 2000, Yi et al. 2005) as well as ob- 
servations that distant early-type galaxies sometimes exhibit 
[Oil] emission lines that may indicate small (r^ 1 Moyr"' or 
less), but significant amounts of star formation (Schade et al. 
1999; WilHs et al. 2002; van Dokkum & ElHs 2003; GSS8; 
Treu et al. 2005). Since cluster early-type galaxies appear 
"red and dead", with little evidence for any residual star for- 
mation, the above explanation would not apply. If verified, 
the small change in cluster bulge colors favors, instead, stel- 
lar population solutions, such as a truncated IMF (Broadhurst 
& Bouwens 2000). 

In summary, our key result is that distant luminous pbulges 
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at redshift z ^ 0.83 have unusually red colors and that such 
red colors appear to be independent of environment (cluster, 
field, local density), pbulge or galaxy luminosity, pdisk color, 
or pbulge fraction (pB/T). Our experiment uses only 21 field 
and 71 cluster galaxies at one redshift and uses one color. This 
initial foray and its conclusions need to be checked by ex- 
panding the sample sizes, broadening the range of colors, and 
deepening the samples. The new surveys should aim to cover 
lower and higher redshifts for both field and cluster galaxies; 
reach lower luminosity bulges especially within later-type spi- 
rals; and achieve sufficient photometric precision to positively 
detect the intrinsic color dispersions. These experiments are 
all now possible with the new, rich, and plentiful HST ACS 
data in the GOODS fields (Giavalisco et al. 2004), Hubble 
Ultra Deep Field (UDF: Beckwith et al. , in preparation), and 
many distant rich clusters (e.g., by the ACS GTO team). With 
the advent of laser guide star adaptive optics on 8-lOm class 



telescopes, stellar population studies of small bulges may now 
also be extended into the near-infrared (thus gaining better 
discrimination among age, metallicity, and extinction) with 
the same spatial resolution of 0.05 arcsec as with HST in the 
optical (e.g., Melbourne et al. 2005). 
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TABLE 1 
Color-Magnitude Fits 



Sample 


No. 


Color" 


Dispersion* 






U-B 


mag 


Cluster - All 


71 


0.438 ± 0.012 


0.10(0.11) 


Cluster - Inner 


23 


0.462 ± 0.021 


0.10(0.12) 


Cluster - Middle 


25 


0.421 ± 0.019 


0.09(0.11) 


Cluster - Outer 


23 


0.436 ± 0.022 


0.10(0.11) 


Field 


21 


0.488 ± 0.029 


0.13 (0.17) 



Note. — 

"Colors are from fits assuming a slope of 0, i.e., that colors are independent of luminosity, using the method of IBeers. Flvnn. & GebhrnHtl fiMI) : errors are 
RMS. ''Dispersions are as measured; values in parentheses are the expected dispersions based on 68% confidence limits derived from GIM2D measurements. 
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Fig. 1 . — Restframe U —B color vs Mb for GSS field and MS 1054-03 cluster pbulges. The solid line aXU -B = 0.44 is set at the median color of the full cluster 
sample. To show any dependence on the environment, the cluster sample is divided into thirds by cluster-centric distance. The symbols show the pB/T ratio and 
their colors indicate pdisk colors. The outlier in the field (boxed) has peculiar morphology. The top-left enclosure shows 68% confidence-limit errors in color 
and luminosity from GIM2D at three luminosities. 



